Quantum correlations in nanostructured two-impurity Kondo systems 
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We study the ground-state entanglement properties of nanostructured Kondo systems consisting 
of a pair of impurity spins coupled to a background of confined electrons. The competition between 
the RKKY-like coupling and the Kondo effect determines the development of quantum correlations 
between the different parts of the system. A key element is the electronic filling due to confine- 
ment. An even electronic filling leads to results similar to those found previously for extended 
systems, where the properties of the reduced impurity-spin subsystem are uniquely determined 
by the spin correlation function defining a one-dimensional phase space. An odd filling, instead, 
breaks spin-rotation symmetry unfolding a two-dimensional phase space showing rich entanglement 
characteristics as, e.g., the requirement of a larger amount of entanglement for the development of 
non-local correlations between impurity spins. We illustrate these results by numerical simulations 
of elliptic quantum corrals with magnetic impurities at the foci as a case study. 

PACS numbers: 03.67.Mn, 73.20.At, 73.20.Hb, 75.30.Hx 
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I. INTRODUCTION 

Along the last decades, the emergence of quantum 
information^ raised a singular outlook in the physical 
science towards the development of a solid-state quan- 
tum engineering 3 This was supported by the contem- 
porary progress of mesoscopic physics, where the mod- 
ern laboratory techniques gave a comprehensive access 
to quantum phenomena in solid-state structures at the 
nanoscale. Presently, many efforts are devoted to the 
design of nanostructured devices for the controlled ma- 
nipulation of light and matter in the search for func- 
tional quantum phenomena. A paradigmatic example is 
quantum entanglement, non-local correlations considered 
a precious resource for quantum-information processing. 

Among the artificial nanostructures stand out the so- 
called quantum corrals. These are arrays of adatoms de- 
posited on the surface of nobel metals, forming a nearly 
closed structure within which surface electrons can re- 
main confined. Quantum corrals have allowed the obser- 
vation of outstanding effects such as quantum mirages, 3 ' 4 
where the introduction of an impurity on one of the 
foci of an elliptic structure forms a ghost image in the 
empty one£ Several authors have analyzed these exper- 
iments with one or more impurities considering different 
configurations.— ~— It was suggested that two impurities 
which are located at the foci of the system will inter- 
act strongly as a consequence of the focalizing properties 
of quantum elliptic corralsj^^— (Alternative locations 
apart from the foci were also suggested^ for the produc- 
tion of quantum mirages). This depends on the com- 
petition between RKKY and Kondo interactions. 19 ' 20 A 
preliminary study of quantum entanglement in such sys- 
tems was done in Ref. [l7L There, it was found that the 
localized spins develop a strong mutual entanglement for 
a small superexchange coupling to the surface electronic 
states, whereas they disentangle for a strong coupling 
(when compared to the corresponding Kondo tempera- 



ture). 

These results suggest that quantum corrals are suit- 
able for essaying possible applications in quantum infor- 
mation and spintronics^i For instance, magnetic impu- 
rities added on the foci could be used as quantum units 
of information (qubits) where quantum correlations (en- 
tanglement) mediated by delocalized electronic states are 
favored by focusing and confinement effects.— ~— Even- 
tually, quantum corrals could serve as elementary proto- 
types for quantum-information processors provided that 
logical gates can be implemented by controlled manipu- 
lation of the interaction between the magnetic qubits^ 

Here we study the entanglement developed between 
two magnetic impurities embedded in a confined elec- 
tronic environment. This corresponds to a nanostruc- 
tured Kondo system, where localized magnetic impuri- 
ties are subject to an effective (RKKY-like^ 3 -) exchange 
interaction mediated by electron-in-a-box states coupled 
antifcrromagnetically (AFM) to the impurities. The key 
element in our discussion is the electron confinement. 24 
Here, in contrast to extended systems with electrons or- 
ganized in a Fermi sea^ the strength of the quantum 
correlations is determined by the electron filling. For an 
extended Kondo model, the ground state of the composed 
two-impurity/electron system forms a spin singlet^ Cho 
and McKenzie^ have shown that this reduces the two- 
impurity subsystem to a (mixed) rotationally invariant 
Werner stated This is a particular family of (mainly 
mixed) states characterized by a single parameter identi- 
fied with the impurity-spin correlation function, defining 
a one-dimensional phase space running continuously from 
a fully uncorrelated classical ensamble to a fully entan- 
gled (singlet) state. We shall see that confined Kondo 
systems with even electronic filling share similar prop- 
erties. An odd electronic filling, instead, takes the two- 
impurity subsystem away from the Werner-state family 
to a region of the Hilbert space with broken rotational 
symmetry and richer entanglement characteristics. Such 
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an extended family cannot be fully described by a single 
observable. Instead, as we shall see, it unfolds a two- 
dimensional phase space. 

Moreover, it is known that entanglement does not 
guarantee non-locality: while the former refers to the 
(non)separability of a quantum state, the latter is defined 
only by the violation of Bell-like inequalities (meaning 
that quantum correlations cannot be modeled by a hid- 
den classical- variable theory). Within the Werner-state 
family, indeed, there exists a subclass of entangled states 
which do not violate Bell's inequalities. The development 
of non-locality requires to overcome a minimal amount 
of entanglement (measured in terms of, e.g., the concur- 
rence^-). We shall see that an odd electronic filling gen- 
erally constrains the impurity-spin system to develop a 
larger amount of entanglement for the violation of Bell's 
inequalities when compared to the case of even electronic 
filling. 

The work is organized as follows. In Sec. [TT] we dis- 
cuss the role played by the electronic filling on the re- 
duced impurity-spin system based on symmetry argu- 
ments, only. In Sec. Mil we introduce a phase-space 
diagram classifying the entanglement properties of the 
impurity-spin system and discuss the conditions neces- 
sary for the development of non-local correlations. Fur- 
thermore, in Sec. HVI we study the complementary entan- 
glement arising between the impurity-spins and a con- 
fined electronic bath. Finally, in Sec. [V] we present nu- 
merical results corresponding to Kondo-model simula- 
tions for an elliptic quantum corral discussing the dis- 
tribution of solutions over the entanglement phase-space 
diagram. 

II. THE INFLUENCE OF THE ELECTRONIC 
STATES ON THE IMPURITY SPIN SYSTEM 

A. Full (pure) electron/impurity state 

We start by considering a generic pure state of the 
composed system ABC with a definite ^-projection of the 
total spin S^, where A and B name the impurity spins 
and C identifies a set of N electrons confined within an 
arbitrary nanoscopic system (see Fig. [1] for the particular 
case of a quantum corral): 

|*) = tt>l*i>+oa| n)l$2)+a 3 | ;t)l$3)+a4| W)l*4>. 

Here, the first ket of each term represents the two- 
impurity spin subsystem while the |$„) are TV-particle 
electronic states of the form 

l*l>=E^nl7n)i (2) 
n 

l$2>=$>M, (3) 

n 

l*3>=X>»frn>. (4) 
n 




FIG. 1: Nanostructured Kondo system: The figure depicts 
two magnetic impurities A and B on the foci of an elliptic 
quantum corral. Confined electronic states C of given filling 
couple antiferromagnetically to the impurity spins, develop- 
ing either local Kondo states or RKKY-like superexchange 
between impurity-spins depending on the coupling strength 
J. 



l*4>=5>£l7;J->- (5) 

n 

The sets {|7^)}, {|7°}}> and {|7^)} are orthonormal 
bases of ./V-electron states with z-projection of the spin 
= — 1, Sp, and Sf + 1, respectively, where Sp 
(total z-component of the spin, including impurities) has 
a definite value. The electronic states |$„) are such that 
($i|$n) = Sin and ($ 4 |$n) = hn, while ($ 2 |$3> + in 
general since both states belong to the same spin sub- 
space. Without loss of generality, in Eq. ([T]) we assume 
real and positive probability amplitudes a n such that 

ErX = L 

The symmetry properties of upon spin rotation 
depend on the electronic filling N, a constraint for the 
total spin. For an even electron number, the ground 
state of the antiferromagnetic Kondo model for the com- 
posed ABC system is a spin singlet, invariant under joint 
spin rotation i 25 i 26 The situation changes for odd elec- 
tron number, where total spin is a half-integer. This is 
reflected in the reduced two-impurity spin system AB 
with striking consequences on quantum correlations, as 
we shall discuss. 



B. Reduced two-impurity spin density matrix 

The reduced (mixed) 4x4 density matrix for the local- 
ized spins pab is obtained by tracing out the electronic 
degrees of freedom C from the full (pure) density matrix 
p = |*)(*|. In the ordered basis {| ft), I U), I it), I W)} 
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that yields 



PAB = 




We first notice that pab does not mix impurity states 
with different spin projection. This results from assum- 
ing that both the electron number and the projection 
of the total spin along the z axis are good quantum 
numbers. Further consequences follow from considering 
those symmetry constraints present in (though not lim- 
ited to) elliptic quantum corrals with magnetic impuri- 
ties at the foci. The reflection symmetry along the minor 
axis implies that pab must be invariant under magnetic 
impurity interchange. This means pab = Pba, with 
a 2 = a 3 and (<S> 3 |<E>2) = (^3^2)* = cos vj, where the 
phase < ip < ir parametrizes the overlap between |$2) 
and l^s). Under such symmetry constraints, Eq. (j6|) re- 
duces to 

/ 1 - 2a\ - a\ \ 

PAB - a 2 2 cosif a\ ' (7j 

V a\ J 

with < a? = 1 - 2af - a\ < 1, < a| < 1/2 and 
< o| < 1. Notice that orthogonal electronic compo- 
nents |$2) and |$3) (cosy = 0) would lead to vanishing 
coherence (off-diagonal) terms in pab- This would pre- 
vent the impurity spins from developing mutual quantum 
correlations such as entanglement (see Sec. Mil for a de- 
tailed discussion). 



By observing the central 4x4 block of matrix ([7]), we 
notice that pab can be expanded as 

Pab = (l-2a 2 2 -al)\n)(n\ + al\U)(U\ (8) 
+ y ^ (l + ^cos^l^X^I, 

5=±1 

where |i/' ± ) = (| t4-)±l It) )/v / 2 are fully entangled triplet 
(+) and singlet (— ) Bell states. Interestingly, pure triplet 
and singlet states are possible for a\ = 1/2 (implying 
a\ = 04 = 0) provided that cost/? = ±1, respectively. 



C. Hilbert-Schmidt decomposition and 
Werner-state component 

It is helpful to rewrite Eq. ([7]) by introducing the 
Hilbert-Schmidt decomposition^ 

Pab = 1 Ta P a A ® a B ( 9 ) 

a,j3— 0,x,y,z 

with physical observables 

r a0 = Tr(a%a B p A B) G R, (10) 
where cr? = tj and erj are the 2x2 identity and Pauli 
matrices, respectively (j = A,B and v = x, y, z). By as- 
suming invariance under magnetic impurity interchange 
and spin rotation symmetry around the z axis, Eq. @ 
reduces to 



l+r z 



2r 02 



PAB = i 





1 V 'zz ^Txx 









r zz - 2r 02 



(11) 



where we find r xx — r yy — 2a\ cosy, r zz = 1 — 4<22, and 
tq z = 1 — 2(a| + a|) by direct comparison with Eq. (UJ. 

When the interaction between electron and impurity 
spins is antiferromagnetic and isotropic, the ground state 
of the composed ABC system, I^g), is either a spin sin- 
glet or a doublet according to the electronic filling (even 
or odd, respectively). In any of these cases the reduced 
density matrix of Eq. (JTTJ) satisfies r xx = r zz = (an 
r X x 7^ r zz , instead, is still possible for excited states). 
Under such conditions we can conveniently write 



Pab = Pab + a ab (12) 



with 

Pab = \ (l 4 x4 + r ]T a%®o- B ) , (13) 

\ a—x,y,z / 

Aab = r -^-{<j A ®t B + t A ®<y Z B). (14) 

Here, p^ B is a Werner-state density matrixi 25 ' 27 This is 
most commonly written as 

PIB = ^14X4 + ^l|*-><*-|, (15) 

parametrized by the so-called singlet fidelity p s = 
(*~|/>Ab|# - ) = (l-3r)/4with0 < p s < 1. Remarkably, 
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p AB is determined by one single observable: the spin-spin 
correlation function (Sa • Sb) = 3r/4 = 1/4 — p B , with 
-3/4 < (Sa • Sb) < 1/4. The A A b of Eq. in- 
stead, is a traceless contribution to p^B (namely, it is 
not a density matrix itself) determined only by the z- 
projection of the impurity spin, since tq z — (S\ + S B ). 
The pab is then fully determined by two observables, 
only: (Sa ■ Sb) and (S A + S B ). This is a consequence of 
the setting r xx = r zz . Otherwise, (S^-Ss) and (S A + S B ) 
would not be sufficient to determine pab completely, pro- 
viding only a partial tomography of the quantum state. 
For (S A + S B ) = 0, pab reduces to a Werner state. This 
is the case when the full (impurity/electron) state is a 
singlet, as discussed in Ref. 25. A (S A + S B ) ^ takes 
Pab away from the Werner-state subspace modifying the 
entanglement features (see Sec.QlT]for a detailed discus- 
sion). 

The setting r xx — r zz also introduces some constraints 
on the electronic parameters of Eq. ([7]). By noticing that 
r xx = 2a\ cos ip and r zz = 1 — ^a\, we find cos ip = —2 + 
l/2d2 with 1/6 < a\ < 1/2. This allows the development 
of pure singlets \ip~) for a\ = 1/2 at the same time it 
forbids pure triplets \ip+) (see Eq. ©). For a\ = 1/4, 
instead, all coherence terms in Eq. ([7]) vanish together 
with cos ip. 

D. Relation between (Sa • Sb) and {S% + S%) 

While the observables (Sa • Sb) and (S A + S Z B ) are 
in principle independent, the introduction of symmetry 
constraints can restrict such freedom. However, symme- 
try alone is not sufficient for determining the ultimate 
relation between them. This, as expected, is eventually 
given by the full Hamiltonian. 

From Eqs. © and (HD we find (Sa • S B ) = (2r xx + 



r zz )/4 = 1/4 + (cosp - l)a| and (S% + S Z B ) = r 0z = 
1 — 2(a\ + a|). The condition r xx — r zz yields 

(S'a + SI) = l + l(S A -S B )-2at (16) 
= -l-l(S A -S B )+2al. 

This equation is obtained from the application of sym- 
metry constraints, only. From it we can identify sev- 
eral scenarios. For the symmetric case a\ — a\ we find 
(S% + S Z B )=0 while (S A ■ Sb) runs freely from -3/4 to 
1/4. This corres pon ds to the Werner state pab = Pab 
discussed in Ref. \2M (see Eq. (fT2|l). We further see that 
when either a\ or a\ equals zero the relation between 
(S^-Sb) and (S A + S B ) must be linear. Otherwise, they 
are free to develop a more general non-linear behavior 
which must be determined by solving the Hamiltonian 
equation. We illustrate this in the following section. 
III. ENTANGLEMENT IN THE REDUCED 
SPIN SYSTEM 

A. Concurrence and phase-space representation 

A standard measure of the entanglement between two 
1/2 spins or qubits is the concurrence C<2& This runs 
from zero for separable (disentangled) states to one for 
maximally entangled (Bell) states. For the reduced two- 
impurity spin density matrix p AB of Eq. Ijlip we find a 
concurrence 



C(pab) = 




(17) 

This can be fully determined by (Sa • Sb) and (S A + S B ) 
provided that r xx — r zz , reducing to 



n/MB> = n,ax<!KS A -S B >|- W(|< S ^ • S *> + i) -MS A + S B y, 0}. CUM 



In Fig. [2] we depict a phase-space diagram classifying 
the entanglement properties of the pab given in Eq. (|12[) 
as a function of (Sa • Sb) and (S A + S B ). The shaded 
area corresponds to the entangled phase with C > 0, 
according to Eq. (TT5)) . This region is defined between 
curve 2: (S A -S B ) = -1/4+ (l -y/1- 3(S A + S Z B ) 2 ) /2, 
and curve 3: (S A ■ S B ) = -3/4 + (3/2)\(S% + S Z B )\, 
along which C = and C = \2\(S Z A + S B )\ - 1|, respec- 
tively. The region defined above curve 2: (S^ • Sb) > 
-1/4 + (1 - ^/1-3(S% + S Z B ) 2 ) /2, encloses the sub- 
space of separable states (C = 0). The region be- 
low curve 3: (S A ■ Sb) < -3/4 + (3/2)\(S z A + S Z B )\, 



instead, corresponds to unphysical states with imagi- 
nary concurrence.— For (S% + S B ) = 0, Eq. (JTHJ re- 
duces to the well-known concurrence for Werner states 
C(pJ B ) = max{2p s - 1 = -(1/2 + 2(8^ ■ S B », 0} (see 
Sec. Ill C[) . This family of states is represented in the 
phase-space diagram of Fig. [5] by curve 1 (dashed line) , 
where the value (Sa ■ S B ) C = —1/4 (p c s = 1/2) defines a 
critical point separating entangled from product Werner 
states. 

Maximum entanglement (C = 1) is reached only at a 
single point ({S% + S%) = 0, (S A ■ S B ) = -3/4) cor- 
responding to the singlet \ip~), from which C decreases 
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-1 -2/3 -1/2 1/2 2/3 1 



(S~+S' B ) 

FIG. 2: Phase-space diagram for the reduced impurity-spin 
system AB: The figure classifies the entanglement properties 
of pab according to Eq. (|18p . parametrized by the observables 
(S% + S Z B ) and (Sa ■ Sb). Curve 1 along (S A + S%) = cor- 
responds to Werner states. Curve 2 defines the critical limit 
separating entangled from disentangled states, along which 
the concurrence C vanishes. Curve 3 marks the limit for 
physically possible states compatible with the symmetry con- 
straints. Curve 4 indicates a subspace of entangled, mixed 
triplet states. Curve H corresponds to numerical simulations 
on an elliptic quantum corral for an odd filling. 



monotonously (eventually vanishing) as (S A • Sb) ap- 
proaches zero. From there on, C increases again up to 
C = 1/3 - ^4/9 - (S Z A + S Z B ) 2 for (Sa • Sb) - 1/4 and 
1/V3 < \(S% + S%)\ < 2/3 (segment 4 in Fig.©, corre- 
sponding to mixed triplet states (i.e., incoherent super- 
positions of |V> + )>| tt)j an d | 11))- We then conclude 
that, remarkably, a (S A + S g ) ^ allows for entangle- 
ment beyond the limit imposed to Werner states, namely, 
C > for (S A ■ S B ) > -1/4. 



B. Violation of Bell inequalities 

The presence of quantum entanglement does not guar- 
antee the development of the non-local quantum cor- 
relations present in maximally entangled (Bell) states 
as, e.g., the spin singlet. Non-locality is prop- 
erly identified through the (violation of) Bell-CHSH 
(Clauser-Horne-Shimony-Holt) inequalities. 31 Some en- 
tangled (non-separable) states can certainly satisfy Bell- 
CHSH inequalities: a well-known example are the Werner 
states with -3/4\/2 < (S A • S B ) < -1/4 (1/2 < p s < 
(1 + 3/%/2)/4). 32 This means that, in such states, the 
statistical properties ascribed to the entanglement can 
be reproduced by using a hidden variable model, i.e., 
quantum features can be modeled classically. 

Following the criteria introduced in Ref. l32l we find 
that the pab of Eq. (IT21 violates the corresponding 
Bell-CHSH inequality for (S A ■ Sb) < -3/4-\/i Sur- 
prisingly, this upper bound coincides with the result 




-1 -1/2 -0.146 0.146 1/2 



<s;+s;> 

FIG. 3: Locality vs. non-locality in the reduced impurity- 
spin system: The value (Sa ■ Sb) = — 3/4V5 ~ —0.53 deter- 
mines a critical limit separating entangled states subject to 
non-local correlations (violating Bell-like inequalities) from 
those compatible with a local hidden- variable model (satisfy- 
ing Bell-like inequalities). The minimal amount of entangle- 
ment required for the development of non-local correlations 
increases with | (S A + S%)\. 

found for Werner states. However, the presence of a fi- 
nite (S A + S Z B ) has some consequences on the minimal 
amount of entanglement (measured in terms of concur- 
rence) needed for violating the inequality. By setting 
(S A -S B ) = -3/4^2 « -0.53 in Eq. (IHJ), we find 

C Ws B >=^ = ^ - \^-^)-^ + S ' b)2 - 

(19) 

Here we see that the minimal concurrence required for vi- 
olating the Bell-CHSH inequality runs fromC = 3/2V2- 
1/2 - 0.561 for Werner states {(S% + S Z B ) = 0) to 
C = 1/V2W 0.707 (| (S Z A + S%)\ = (l-l/V2)/2 « 0.146), 
see Fig.[3l This means that states with a finite (S A + S B ) 
require a larger amount of entanglement for violating lo- 
cality than Werner states. 



IV. ENTANGLEMENT BETWEEN IMPURITY 
SPINS AND CONFINED ELECTRON 

Entanglement in a quantum system depends on the 
partition choice, namely, on the particular (pair of) sub- 
systems under study. We now explore the entanglement 
developed between the impurity-spins and the electronic 
subsystems. A good measure of this entanglement is the 
von Neumann entropy E of the reduced impurity-spin 
density matrix p AB ^^ defined as 

4 

E(pab) = - ^ w " log 2 ( 20 ) 

n=l 

where 

wi = l + l(S A -S B ) + l(S z A + S B ), 



6 



^ = - + -(S A -S B )--(S Z A + S Z B ), 

0J3 = \ -1(Sa-S b ) + 1\(Sa-S b )\, 
11 2 

w 4 = |-g<SA-Ss)--|(SA-Ss)|, 

are the eigenvalues of pab- The E(pab) runs from 
zero for separable states to two for maximally entan- 
gled states. Absolute maxima and minima of E(pab) 
are found along the curve defined by (S A + S B ) = 
at points (S^ • S B ) = and (S A ■ Sb) = —3/4, re- 
spectively (see Fig. [4j. This coincides with what was 
found in the case of global impurity-spin/electron sin- 
glet states^ where E(pab) = 2 indicates the formation 
of local Kondo singlets between impurity spins and elec- 
trons. E(pab) = corresponds to the development of an 
impurity-spin singlet separated from the electronic sys- 
tem. The introduction of a finite (S A + S Z B ) produces a 
decorrelation between impurity-spins and electrons, re- 
ducing the amount of entanglement E(pab) as shown in 
Fig. 21 Such decorrelation, however, is not complete (ex- 
cept in the vicinity of the impurity-spin singlet), so that 
we can generally claim that the impurity-spin subsystem 
shows a relatively high correlation with the electronic de- 
grees of freedom. 

Complementary information can be obtained by quan- 
tifying the entanglement of one single impurity-spin with 
the rest of the system. This can be done by calculat- 
ing the von Neumann entropy E of the reduced density 
matrix of the impurity spin pa (equivalent results are 
obtained for the impurity spin B). We find 

EiPA) = - l±M) log2 l±M) (21) 

1-(S% + S B ) 1-(S% + S B ) 
2 l0g2 2 ' 

where pa = {&a + (&A + ^b)^)/^ 1S independent of the 
correlation function (S^ • S^). The E(pa) is bounded 
between zero (uncorrelated spin-impurity) and one (fully 
entangled spin-impurity). Full entanglement of single 
impurity-spins is reached for (S A + S B ) = 0, in agreement 
with previous results for global impurity-spin/electron 
singlet states^ The presence of a finite (S A + S B ) decor- 
relates the single impurity spins from the rest of the 
system. Full decorrelation (separation) would require 
(S A + S B ) = ±1. However, symmetry constraints impose 
(S Z A + S Z B ) < 2/3 (see Fig. [2]) forcing the entropy E{p A ) 
to be lower-bounded by a finite value close to 0.65. This 
means that single impurity spins remain highly correlated 
with the rest of the system even for odd electronic filling. 

V. NUMERICAL SIMULATIONS ON 
QUANTUM CORRALS 

The results discussed above are valid for general two- 
impurity models interacting with a confined electronic 
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FIG. 4: Entanglement between confined electrons and im- 
purity spins given by the von Neumann entropy E(pab) of 
Eq. Q20p parametrized by the observables (SJ + Sg) and 
(Sa ■ Sb). Broken- line and solid curves correspond to nu- 
merical simulations on an elliptic quantum corral for even 
and odd electrons respectively. 

bath. We now perform a case study of the entangle- 
ment features for two impurities located at the foci of 
an elliptic corral confining surface electrons. Each impu- 
rity interacts antiferromagnetically with the spin of the 
band electron located at the corresponding focus (Fig. 
[1}, modeled by the Hamiltonian: 

H = H C + J(S a .<ta + S b .ct b ), (22) 

with J > 0, 

St.at = Slot + \{St.ai + Sr.af), (23) 

af = c^Cii, and of = (n it - n^)/2. Here, n ia and c ia 
are, respectivelly, the number and destruction operators 
for electrons with spin a on focus i = A,B. In the basis 
of eigenstates |a) of the ellipse, these local operators can 
be expanded as a a = J2 a ^aiC a <j where c aa and *& a i are 
the destruction operator and amplitude in state \a). In 
this basis the electronic spin operators are expressed as: 

4 = J2 *ai<*«a<4 I t<Wl.- ( 24 ) 

The H c in Eq. (|22|) accounts for electron confinement in 
an ellipse with hard walls (eccentricity e = 0.6 in numer- 
ical simulations)^ 

We perform many-body simulations of Eq. (|22j) ap- 
plying numerical techniques as exact diagonalization or 
Lanczos, considering up to 10 electronic levels in the el- 
lipse and different particle fillings N (either even or odd) . 
For odd electronic fillings we consider just the = 1/2 
subspace. Results for S% = —1/2 are equivalent due 
to up-down spin symmetry (something reflected by the 
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left-right symmetry of Figs. [2 [3] and [4]). We also chose 
the Fermi energy to coincide with the 23 rd eigenstate 
1^23), reproducing the experimental conditions of Ref. H 
This configuration is not crucial for the main conclusions: 
other eccentricities and Fermi levels show similar features 
as long as the wave function bears an appreciable weight 
at the foci of the ellipse. 

We first follow Eq. ([To]) by calculating numerically the 
spin correlation function, (Sa ■ Sb), and the z-projection 
of the localized spins, (S A + S Z B ), for different values of 
the coupling strength J. We find a coefficient 04 = 
for small odd fillings (N < 3) and a 04 ^ otherwise, 
leading to a non-linear relation between the observables. 
Results for 10 levels and 9 particles are plotted in Fig. 
[2] (curve H represented by a broken line). By following 
its route along the phase-space diagram one observes the 
entanglement features of the impurity spin subsystem un- 
der such conditions, running from separable to entangled 
areas within the (Sa ■ Sb) > region. More complex 
behaviors are to be expected for less restrictive symme- 
tries as shown in Ref. HH- Results for even fillings are 
indicated with a broken line coinciding with curve 1. 

Finally, we quantify the entanglement developed be- 
tween the localized impurities and the electronic bath in 
the ellipse based on Eq. (I2TJ1) . The results are depicted in 
Fig. HI the broken line for an even filling (10 levels and 
10 particles) and the full line for an odd filling (10 levels 
and 9 particles). Here we notice that the entanglement 
monogamy^ is clearly satisfied by observing that a large 
entanglement entropy between impurity spins and elec- 
trons corresponds to a low concurrence between impurity 
spins, and vice versa. 

VI. CONCLUSIONS 

We have identified the ground-state entanglement 
properties of two magnetic impurities (1/2 spins or 
qubits) interacting with confined electronic backgrounds 
by composing a phase diagram for the classification of 
separable, entangled, and unphysical states. This was 
done by studying the impurity-spin concurrence as a 
function of the spin correlation function and the z- 
projection of the localized spins. Such quantum cor- 



relations show a striking dependence on the electronic 
filling: while an even filling reduces the impurity spin 
system to a highly symmetric Werner state, an odd fill- 
ing breaking spin rotation invariance give rise to richer 
features. We also discussed the relation between entan- 
glement and (non)locality by studying the conditions for 
the violation of Bell inequalities in the impurity spin sys- 
tem. Additionally, we obtained complementary informa- 
tion by studying the entanglement developed between the 
the impurity spins and the electronic background. 

As a case study, we performed corresponding numeri- 
cal simulations for elliptic quantum corrals consisting of 
a pair of impurity spins located at the foci coupled an- 
tifcrromagnetically to a background of surface electrons 
confined by the elliptic potential. By increasing the in- 
teraction parameter, there exists a transition from an 
RKKY-like regime (where the impurity spins are mutu- 
ally entangled) to a localized Kondo-like regime (where 
impurity spins decorrelate, developing a strong entangle- 
ment with the electronic environment). We notice that 
electron confinement may help to preserve quantum cor- 
relations by protecting the impurity spin systems from 
external fluctuations. Moreover, impurity spin correla- 
tions could be accessible experimentally by implement- 
ing magnetic STM tips or similar devices. We also ex- 
pect these results to be valid for other kinds of back- 
grounds such as spin chains. Quantum corrals offer an 
ideal scenario to study entanglement between impurities 
and corral wave functions or other possible partitions. 
This quantum information interpretation is an alterna- 
tive way of analyzing crossovers between different quan- 
tum states and sheds light onto other relevant properties 
stemming from focalization and confinement. 
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